Chiral Magnetic Wave at finite baryon density and 
the electric quadrupole moment of quark-gluon plasma in heavy ion collisions 
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Chiral Magnetic Wave (CMW) is a gapless collective excitation of quark-gluon plasma in the 
presence of external magnetic field that stems from the interplay of Chiral Magnetic (CME) and 
Chiral Separation Effects (CSE); it is composed by the waves of the electric and chiral charge 
densities coupled by the axial anomaly. We consider CMW at finite baryon density and find that it 
induces the electric quadrupole moment of the quark-gluon plasma produced in heavy ion collisions: 
the "poles" of the produced fireball (pointing outside of the reaction plane) acquire additional 
positive electric charge, and the "equator" acquires additional negative charge. We point out that 
this electric quadrupole deformation lifts the degeneracy between the elliptic fiows of positive and 
negative pions leading to V2{'k^) < V2{tt~), and estimate the magnitude of the effect. 

PACS numbers: 11.40.Ha,12.38.Mh,25.75.Ag 
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Introduction. — The axial anomaly has been found to 
induce the following two phenomena in the quark-gluon 
plasma subjected to an external magnetic field: the Chi- 
ral Magnetic Effect (CME) and the Chiral Separation 
Effect (CSE). The CME is the phenomenon of electric 
charge separation along the axis of the applied magnetic 
field in the presence of fluctuating topological charge 
[ll-Q. The CME in QCD coupled to electromagnetism 
assumes a chirality asymmetry between left- and right- 
handed quarks, parametrized by an axial chemical poten- 
tial fiA- At finite iia, an external magnetic field induces 
the vector current ji — Tp^ttp: 
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in our present convention the electric current is ejv ■ 

Recently, STAR % 0] and PHENIX 0, Q Collabora- 
tions at Relativistic Heavy Ion Collider at BNL reported 
experimental observation of charge asymmetry fiuctua- 
tions possibly providing an evidence for CME; this inter- 
pretation is still under intense discussion, see e.g. 13, Hi 
and references therein. 

The Chiral Separation Effect (CSE) refers to the sep- 
aration of chiral charge along the axis of external mag- 
netic field at finite density of vector charge (e.g. at finite 
baryon number density) [l3l | . 



The resulting axial cur- 



rent is given by 
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where /iy is the vector chemical potential. Both CME 
and CSE effects have been proved robust in holographic 
QCD models in strong coupling regime_llj-|20| as well 
as in lattice QCD computations 2l|, |22|. The effects 
also persist in relativistic hydrodynamics, as shown in 
Ref.fH. 



Recently, two of us studied the properties of the " Chi- 
ral Magnetic Wave" (CMW) |2J] stemming from the cou- 
pling of the density waves of electric and chiral charge 
induced by the axial anomaly in the presence of external 
magnetic field; a related idea has been also discussed in 
[isf . The CMW is a gapless collective excitation; its ex- 
istence is a straightforward consequence of the relations 
eq.iH]) and eq.(I2]). Indeed, consider a local fluctuation of 
electric charge density; according to eq.(I2|) it induces a 
local fluctuation of axial current. This fluctuation of ax- 
ial current in turn induces a local fluctuation of the axial 
chemical potential, and thus according to eq.([T]) a fluc- 
tuation of electric current. The resulting fluctuation of 
electric charge density completes the cycle leading to the 
CMW that combines the density waves of electric and 
chiral charges. 

The plasma created in heavy ion collisions possesses a 
finite baryon density. The CSE [l2, 13, 23j then implies 
the separation of chiral charge: the "poles" of the fireball 
acquire the chiral charges of opposite sign. The CME 
current at the opposite poles then according to eq.(|Tl) 
fiows in opposite directions, as argued recently in [26|. 
In this letter we will show that CMW induces a static 
quadrupole moment of the electric charge density. 

Chiral Magnetic Wave. — The Chiral Magnetic Wave 
(CMW) is a long wavelength hydrodynamic mode of chi- 
ral charge densities; their propagation in space-time is 
described by the following equation [24] 



(3) 



where is the wave velocity and Dl is the longitudinal 
diffusion constant. 

In the case of Nf quark flavors with electric charges g/ 
there will be Nf independent CMWs with the velocities 
and longitudinal diffusion constants determined by g/, 
eB and T. In this paper we consider the propagation of 
u and d flavored CMWs, since there is no net density of 
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strange quarks in the plasma. The full flavor symmetry 
U{2)f contains J7(l)„ x U{l)d which defines independent 
U{1) flavor symmetries of u and d quarks. Considering 
the same triangle anomalies leading to CME and CSE 
that now involve each of these U{1) symmetries, one ob- 
tains 



(4) 



where fi^ are chemical potentials of U{l)f. From the 
results of 12J] and (jH) we then derive Nf independent 
CMWs of flavored chiral charge densities with ve- 
locities given by 
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We obtain and from the computation in Ref. 24 1 
performed in the framework of Sakai-Sugimoto model in 
the large Nc quenched approximation. Each quark of fla- 
vor / interacts with the magnetic field of effective mag- 
nitude QfcB, we replace eB with qfcB in the arguments 
of % and Dl as functions of eB: 



v^ieB^qfcB), D{ = Dl {eB ^ qjeB) . (6) 

We neglect the isospin asymmetry between the u and d 
flavors at the time of plasma creation, and take the equal 
initial chemical potentials ^ly = /x^/S. The shape of the 
initial " almond" of QCD matter produced in a heavy ion 
collision is taken by using the phenomenologically suc- 
cessful KLN model [2^ based on parton saturation and 
kx factorization. Au — Au collisions have been simulated, 
with realistic Woods-Saxon nuclear densities. The axial 
chemical potentials at the initial time are set to zero. 

We then solve the CMW equation numerically and 
find that it generates the separation of chiral charge, as 
shown in Fig[T]- the quark-gluon plasma acquires a "chi- 
ral dipole moment" . 

We evaluate the total electric charge distribution by 
super-imposing the waves of different flavors weighted by 
their charges. 
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The resulting distribution is shown in Figl2] for clarity, 
we have subtracted the charge density distribution with- 
out the CMW. As argued above qualitatively, the quark- 
gluon plasma indeed acquires an electric quadrupole mo- 
ment. The "poles" of the produced fireball (pointing 
outside of the reaction plane) acquire additional pos- 
itive electric charge, and the "equator" acquires addi- 
tional negative charge. It is very important to note that 
this pattern of charge separation does not depend on the 
orientation of magnetic field. This means that the effect 
should survive even after the event averaging. 




FIG. 1: Chiral charge density in the plane transverse to the 
beam axis; magnetic field strength eB — m%, lifetime of mag- 
netic field T = IQ fm, temperature T = 165 MeV, impact 
parameter fe = 3 fm. 




FIG. 2: Electric charge density in the transverse plane (back- 
ground subtracted, see text); same parameters as in Fig. 1. 

From the Electric Quadrupole Moment to Charge- 
Dependent Elliptic Flow. — The expansion of the quark- 
gluon plasma produced in heavy ion collisions is char- 
acterized by strong anisotropic collective fiow driven 
by the gradients of pressure that transforms the spa- 
tial anisotropy of produced matter into the momentum 
anisotropy of the produced hadrons. Since the fireball 
of quark-gluon plasma produced in an off-central heavy 
ion collision has an elliptical almond-like shape, the gra- 
dients of pressure make it expand predominantly along 
the minor axis, i.e. in the reaction plane - this is the 
"eUiptic flow" (for a review, see [28|). As a result, the 
electric quadrupole deformation of the plasma described 
above will increase the elliptic flow of negative hadrons, 
and decrease the elliptic fiow of positive hadrons leading 
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< ■ However, the large differences in the ab- 
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sorption cross sections of antiprotons and protons, and 
of negative and positive kaons in hadronic matter at finite 
baryon density are likely to mask or reverse this differ- 
ence in the hadron resonance "afterburner" phase of a 
heavy ion collision. On the other hand, the smaller dif- 
ference in the absorption cross sections of negative and 
positive pions potentially may make it possible to detect 
the electric quadrupole moment of the plasma through 
difference of elliptic flows of pions, V2{tt^) < V2{tt^). 

Let us now quantify this statement and estimate the 
magnitude of the effect. The net electric charge density 
is a sum of the initial charge density and the modulation 
caused by the propagating CMW. Keeping the leading 
term in the multipole expansion, we thus write 



Pe = Pe - 2ge COs(20) = Pe 
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Note that the quadrupole term ^ only re-shuffles the 
initial charges such that the net charge density pe is larger 
out-of-plane than in-plane. 

The bulk evolution of the system should not be sig- 
niflcantly affected by the quadrupole deformation as Qe 
is small compared to the overall density. We therefore 
assume that the overall elliptic flow is not modified and 
the charged hadron distribution in the azimuthal angle 4> 
with respect to the reaction plane is still given by 
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= iVo [1 + 2^2 cos(2<?!))] . 
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Nonzero pe translates into the difference of the yields 
of positive and negative charges. Let us split the to- 
tal charged hadron number at hadronization Nq into 
the numbers of positive and negative hadrons: Nq ~ 
+ N- . With zero pe , one has — iV_ . An isotropic 
Pe stemming from the baryon stopping in a heavy ion 
collision leads to non-zero difference (iV+ — N^) > 
independent of the azimuthal angle cj). When there is, 
in addition, a quadrupole deformation — 2(7e cos(20), the 
asymmetry (N-^ — N^) develops an azimuthal angle de- 
pendence. In the presence of a strong radial fiow, the 
direction of the momentum of a fluid element coincides 
with the direction of the radius-vector pointing to the 
position of this fluid element. Therefore, a quadrupole 
deformation in coordinate space should translate into 
a quadrupole deformation in momentum space, and we 
thus expect a non-zero 0-dependent difference between 
iV+ and N_ : 



N+ = N+- NgC0s(2(j}); 
= iV_ +iV,cos(2(^). 



The difference now becomes 



1 - 



-2cos(20) 



(10) 
(11) 



.(12) 



While the (iV+ — iV_) is due to the charge density pe, 
the (/)-dependent part is due to the quadrupole deforma- 
tion qe- Since the two components of net charge density 
corresponding to the flrst and second terms in Eq. are 
transformed into the final plus/minus difference via the 
same bulk evolution, it is plausible to assume that the 
ratio between the two components - the quadrupole and 
the background parts in Eq. ([T2|) - remains the same as 
in Eq.®: 
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This translates the electric quadrupole deformation of the 
plasma into the charge dependence of the elliptic flow of 
hadrons. Indeed the 0-dependence in Eq. ([T2l) leads to a 
charge-dependent elliptic flow term; 



= N±[l + 2v2 cos(20)] Ri N± 1 + 2v2 cos(20) 



q^^ /iV^ 
^ \flej \N++N. 
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in deriving this relation we assumed V2 << 1 and 
(iV+ - N^) << N±. The elliptic flow thus becomes 
charge dependent: 



V2 T 



A, 



(15) 



where A± = (iV+ — N^)/{N^ + 7V_) is the charge asym- 
metry in the plasma. 

The magnitude of the effect: numerical simulation. — 
As described above, we have computed the evolution 
of the right and left chiral components of the u and d 
quarks according to equation ([3]) (at zero rapidity) in a 
static plasma. For simplicity, we assume the temperature 
to be uniform within the almond. At the boundary of 
the plasma, the chiral symmetry is broken and therefore 
we set v-^ = 0. In the transverse (w.r.t. the magnetic 
field) direction, we assume a diffusion with a diffusion 
constant Dt estimated j2^ as Dt = (27rr)~^ within 
the Sakai-Sugimoto model. The difference in the elliptic 
flows of positive and negative pions is given, within our 
approximation, by eq. (fT5l) . In Figl3]we present the ratio 
r = 2qe/pe of the electric quadrupole moment 



(16) 
(17) 



qe = j RdR dd, cos{24>) {jl{R, 4,) - flB=oiR, 

to the total charge density 

Pe = J RdR d0j°B^o(i?,0). 



In this computation we took the impact parameter de- 
pendence of magnetic field from with the maximal 
value eB\max = rn^. To convert this ratio into the dif- 
ference of the elliptic fiows of positive and negative pions 
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The normalized electric quadrupole moment r, 
'2 T = 165 MeV. 



FIG. 3: 

according to Eq.([T5|) we also have to estimate the electric 
charge asymmetry A± in the quark-gluon plasma that 
varies between and 1. We do this using the baryon 
chemical potential and temperature at freeze-out ex- 
tracted [29,] from the data, and for the energy of ^/s = 11 
GeV estimate A± ~ 0.3; note that at finite baryon den- 
sity the asymmetry in the plasma and at freeze-out may 
differ. The lifetime of magnetic field in the plasma is still 
uncertain: the initial pulse of magnetic field rapidly falls 
of with time [s*] , 'so'] , but the induction in electrically con- 
ducting quark-gluon plasma was estimated to drastically 
extend the lifetime of the magnetic field according to the 
Lenz's rule, perhaps making it last for the entire lifetime 
of the plasma [31|. Choosing for the sake of estimate 
r = 8 fm and the corresponding r ~ 0.04 for mid-central 
collisions in the plot of Fig. [3]and using A± ~ 0.3 we es- 
timate the difference of tt" and tt"^ elliptic flows in AuAu 
collisions at -i/s = 11 GeV as 



A- 



,CMW _ 



V2{n ) - t;2(7r+) ~ r A^ 



0.01 



(18) 



The pion elliptic flow at ^/s — 8.7 GeV in mid-central 



Pb — Pb collisions is about V2 — 0.03 [3^- Therefore, 
the difference between the elliptic flows of positive and 
negative pious may be as big as ~ 30%; this would clearly 
make it observable. 

Uncertainties and outlook. — The main source of un- 
certainty in our computation is the lifetime of magnetic 
field that will need to be evaluated numerically within 
relativistic magnetohydrodynamics (MHD). Our treat- 
ment of expansion and evolution of the plasma has been 
quite crude, and will also need to be refined by an MHD 
computation. Possible backgrounds to the effect consid- 
ered here have to be carefully analyzed, and include the 
Coulomb interaction of produced pious and the difference 
in absorption of negative and positive pious in dense res- 
onance gas. Since the CMW can propagate only in the 
chirally symmetric phase, the electric quadrupole defor- 
mation of the plasma can provide a signature of chiral 
symmetry restoration in heavy ion collisions. 
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